Nonequilibrium thermal entanglement in three-qubit XX model 
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Making use of the master equation and effective Hamiltonian approach, we investigate the steady 
state entanglement in a three-qubit XX model. Both symmetric and nonsymmetric qubit-qubit 
couphngs are considered. The system (the three qubits) is coupled to two bosonic baths at different 
temperatures. We calculate the steady state by the effective Hamiltonian approach and discuss the 
dependence of the steady state entanglement on the temperatures and couplings. The results show 
that for symmetric qubit-qubit couplings, the entanglements between the nearest neighbor are equal, 
independent of the temperatures of the two baths. The maximum of the entanglement arrives at 
Tl = Tr. For nonsymmetric qubit-qubit couplings, however, the situation is totally different. The 
baths at different temperatures would benefit the entanglement and the entanglements between the 
nearest neighbors are no longer equal. By examining the probability distribution of each eigenstate 
in the steady state, we present an explanation for these observations. These results suggest that the 
steady entanglement can be controlled by the temperature of the two baths. 
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When we talk about a real physical system, we should 
take the effects of its environment into account because 
all quantum systems interact unavoidably with their sur- 
roundings. A quantum system that can not isolate from 
its environment is usually referred to open quantum 
system [l|. The dynamics of open quantum systems 
can be described by quantum master equations in the 
Schrodinger picture, or Langevin equation in the Heisen- 
berg picture • The coupling of environment to a system 
definitely changes the properties of the system, such as 
the geometric phase [3|] and entanglement j4l]. 

Entanglement is a quantum resource that has no clas- 
sical counterpart. It was first recognized in 1935, and has 
been widely studied in recent years due to its potential 
applications in quantum information processing_|5| . The 
environment can either induce entanglement |6| or de- 
crease entanglement (in this case environment often lead 
to a death of entanglement, which is usually called entan- 
glement sudden death Q). When a quantum system is 
in contact with a heat reservoir at a fixed temperature, 
the system will relax into a thermal equilibrium state 
p(T)=e~'^" /TT{e~f^") eventually. The entanglement of 
this state is called thermal entanglement, which has been 
extensively studied in the past decade 

In this paper, we shall study a different type of thermal 
entanglement. It depends on temperature but its state is 
not statistically equilibrium. This situation arises when a 
quantum system interacts with two bosonic baths at dif- 
ferent temperatures. The state of the quantum system 
eventually arrived is not a statistically equilibrium but a 
steady state @ . In Ref . [l0| , the heat transport were stud- 
ied for such a system. In Ref. [ll| , the entanglement of a 
two qubit XX chain were studied for both identical and 
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FIG. 1: A schematic representation of our model. The two- 
level systems are connected to two bosonic baths held at dif- 
ferent temperatures, Tl and Tr, respectively. 

different qubits. In this paper, we will study the entan- 
glement of a three qubit XX chain coupled to two baths 
at different temperatures. Two cases, i.e., symmetric and 
nonsymmetric qubit-qubit couplings, are considered. 

Consider a spin chain consisting of three spins with 
XX interaction. The Hamiltonian for this system has 
the form, 

3 

i/=^ |<+Ji(c7+a2-+-7-4)+J2(4a3-+a2-4), (1) 

n=l 

where cr^ and are Pauli operators for the nth 
spin, and Ji and J2 denote coupling constants. 
Spin 1 and 3 interact with two sperate bosonic 
baths at different temperatures and Tr, re- 

spectively (see FigH]). The Hamiltonian for each 

bath j^L.R is given by Z]„ '^nj^ij^nj and 

the interaction between the spin and its bath is 

described by Hsbj^ct^ J2n 9n^bn^j+a~ J2n dn^^l^j = 

^jtM^J^A'+^j^M^i.M = 1>3 for spin), here the oper- 
ator Xj^fj^ is an eigenoperator of the system Hamiltonian 
satisfying [H, X^^]—±ujj_f^X^^, and Fj,^ acts on the bath 



2 



degrees of freedom. We assume the both baths are in the 
uncorrelated thermal equihbrium states. Then the den- 
sity operator of the bath is pBj=e^'^J^^J /Tr(e^^J^^J ). 
In this paper we set fc^^Ti^l for simphcity. Under this 
assumption, the dimension of the transition frequency e, 
couphng constants Ji , J2 and the temperatures Tl , Tr 
are same. In our following discussion, the value of these 
parameters only stand for the ration relation between 
them. The dynamics of the system affected by the two 
baths can be described by the quantum master equation 
0, which is obtained by tracing out the bath variables 
within the Born-Markovian approximation. This equa- 
tion can usually be arranged in Lindblad form 



d_ 

dt 



p = -i[H, p] + ClP + CbP, 



(2) 



where ClP {Crp) is the dissipative term due to the cou- 
pling of the system to the left (right) bath. In order 
to find an exact form of the dissipative term, we go to 
the bases composed by the eigenstates \si) of the system 
Hamiltonian with the corresponding eigenvalues Xi 



|si) = |lll),Ai--e, 
|s2) = |000),A2=-^e, 

|s3) = ^(sin0|llO>-|lOl) + cos0|Oll)),A3=-e-J, 

|s4)=cos6i|110)-sin6i|011), A4=^£, 

l'S5) = ^(sin6i|110) + |101) + cos6i|011)),A5 = -e+J, 

|s6) = ^(cos6i|100)-|010) + sin6i|001)),A6=--e-J, 
|s7)=sin6i|100)-cos6'|001), A7=-ie, (3) 

|sg) = ^(cos6i|100> + |010> + sin6i|001)),A8=--e-l-J, 



where J— Jf + J2, and tan 6*=^. In this representa- 
tion, the dissipative term Cjp {j^L, R) can be written 
as 
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with eigenfrequcncies, 

LJi = e + J, UJ2 ^ £, CJ3 ^ e - J, 

and corresponding eigcnoperators 

V2 1 
= — cos 6*153) (si I - ■^\s6){si\ 

- — C0S6'|S7)(S5| -I- — C0S6I|S2)(S8|, 



(5) 
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L,2 



S4)(si| - sin6'|s6)(s3| 



-^L,3 — 



+ sin0|s8)(s5| + sin0|s2)(s7|, 
— cos 6*155) (si I + — cos 6*137) (S3 I 
1 \/2 

+ 2l*8)(*4l + — C0S6'|S2)(S6|, 



\f2 1 
Xr,\ = — sin6'|s3)(si| -|s6)(s4| 

-h — sin6'|s7)(s5| -f — sin6'|s2)(s8|, 
\/2 

Xr,2 = — |s4)(si| - COS 6*156) (S3 I 

^ C0S6'|S8)(S5| - C0S6'|S2)(S7|, 

-^fl,3 = — sin 6*155) (si I - — sin 6*157) (S3 1 

1 \/2 
- 2l^8)(*4| + — sin6'|s2)(56|. (6) 

The spectral density in Eq.(|3]) is defined by 



j(j')(o 



(7) 



In this paper, the bath is assumed as an infinite set of 
harmonic oscillators, so the spectral density has the form 
j(-')(w^)=7j(a;^)nj(w^), where nj-(w^)=(e'^j"''-l)"i and 
J^^\—LU^)=e^^'^i^ J^^\uj^). On the same footing as 
the Born-Markovian approximation, we take here a 
Weisskopf- Winger form for the coupling constant as 
7j(a'^)=7j, i.e., the coupling constant is spectrum in- 
dependent. 

The master equation can be solved by the fourth-order 
Runge-Kutta method and the steady state can be reached 
when we set the evolution time long enough . Differ- 
ent from this method, we solve the steady state of Ea.Q 
here with the help of effective Hamiltonian approach [13] ■ 
The main idea of this method can be described as fol- 
lows. By introducing an ancilla, we map the density 
matrix and the master equation to a state vector and 
a Schrodinger-like equation. The solution of the mas- 
ter equation can be obtained by mapping the solution 
of the Schrodinger-like equation back to the density ma- 
trix. Rewrite the master equation in the Lindblad form 

[a as 9.p=-^[H,p]-^J:, (LlLkP+pLlLk-2LkpL{). 

Then the effective Hamiltonian is defined by TLt — 

H-H^-^ LlL,~^Ek Lt^Lt+^ ^tU^ where 
A denotes the operator for the ancilla, which is defined 
by {em\0^\en) = {En\O^Em). Here {|i;„)} and {|e„)} are 
time independent bases for the original system and the 
ancilla. The initial state independent steady state cor- 
responding to the eigenstate of Tip with eigenvalue zero. 
Then we can study the steady state properties for such 
a system. 
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FIG. 2: Steady state entanglement C12 and C13 as functions of 
the bath temperatures Tl and Tr. Here e=l and Ji=J2 = 1. 




FIG. 3: Steady state entanglement C12 and C13 as functions 
of the bath temperatures Tl and Tr. Here e=3 and Ji = J2 = l. 



Choosing the concurrence [lj| as the measure of en- 
tanglement, we first study the steady state concurrence 
in the case of symmetric qubit-qubit couplings {Ji—J2)- 
In Figs (21 and [31 we plot the steady state concurrence 
C12 (entanglement between spins 1 and 2), and C13 (en- 
tanglement between spins 1 and 3) at different temper- 
atures Tl and Tr. In FigO we set e=l and in FigU 
£=3. The coupHng strength are chosen to be Ji=J2 = l. 




FIG. 4: (Color online) Left column, the steady state en- 
tanglement as functions of the mean bath temperature 
Tm^{Tl+Tr)/2 and the temperature difference AT^Tl-Tr 
in non-symmetrical case. Right column shows the correspond- 
ing concurrences change with AT in the case of Tm~0.2 
(blue-solid), Tm=0.4 (green-dashed), Tm=0.6 (red-dotted), 
rM=0.8 (cyan-dash-dot), and Tm~1 (pink-circle). Here e=3, 
Ji=0.5 and J2=2.5. 
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FIG. 5: (Color online) The probability distribution for eigen- 
states \si},i=3, 4, • ■ • , 8 in the steady state as a function of the 
temperature difference AT. In top-left figure we set e=l and 
Ji= J2 = l, while in top-right figure we choose £=3, Ji=0.5 and 
J2—2.5. In both two figures, the mean bath temperature Tm 
is fixed with Tm=0.4. In bottom figure, we plot the probabil- 
ity distribution for Iso) with different mean bath temperature. 
The mean bath temperatures correspond to the right column 
of FiglU 




FIG. 6: Steady state entanglement as functions of the cou- 
pling strength Ji and J2 in the case e=3, Tm=0.5, AT= — 0.4. 

Some features can be observed from the figures. First, 
although the value of C12 and C13 are quite different, 
the tendencies of them affected by T^ and Tn are simi- 
lar. The other interesting feature is that, in both Fig[2] 
and FiglSl C12 and C13 are symmetric about T^—Tji, and 
the peak(maximum) appears when TL=Tfj. This means 
that, Ci2=C23 holds in this case. In other words, the 
entanglements between the nearest neighbors are equal 
although the two temperatures are different. Moreover 
the larger the temperature difference, the smaller the en- 
tanglement. 

In nonsymmetric case {Jiy^J2), the results are quite 
different. We first plot the concurrences between any 
two spins as functions of the mean bath tempera- 
ture Tm—^{Tl+Tji) and the temperature difference 
AT=Tl—Tij in FigHl The parameters in the figure are 
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chosen e=3, Ji=0.5 and J2=2.5. It can be easily seen 
from the figure that the peak does not appear at Ar=0 
(For more clearly, see the right column of FigHJ, i.e., for 
a fixed value of mean bath temperature, the peak ap- 
pears at the points where the temperature difference is 
not zero, suggesting that temperature difference benefits 
the entanglement in nonsymmetric case. If one wants 
to get a larger entanglement when the qubit-qubit cou- 
plings are not equal, a specific temperature difference 
is necessary. The other interesting phenomenon in this 
case is that when the difference between the coupling 
constant Ji and J2 is large enough (for example, in our 
figure Ji=0.5, J2=2.5), the difference between C12 and 
C13 becomes negligible. In other words, C13, the concur- 
rence between the next-to-nearest neighbor qubits, tends 
to the concurrence between the nearest neighbors with 
weak coupling (in our condition, C12), and they all are 
smaller than that in the case with stronger couplings (in 
our condition, C23). 

To understand these features, we calculate the prob- 
ability distribution for the eigenstates of the Hamilto- 
nian H in the steady state, this distribution is a func- 
tion of temperature difference AT as shown in FigjSl 
Although the entanglement does not satisfy the role of 
superposition [1^, the probability distribution are helpful 
to understand the features we found in the entanglement. 
In the top-left figure of Figl^l we plot this distribution for 
the symmetric case, i.e., Ji=J2 = 1. The other param- 
eters are e=l, Tm=0.4. From the figure, we find that 
the probability distributions for all the eigenstates are 
symmetrical about AT=0, and |s6)J'S3) dominate the 
distribution. The distribution for \si) and |s2) are not 
shown in this figure, because these two states are sepa- 
rable. On the other hand, for symmetric couplings we 
have tan 0=1, then the eigenstates can be divided into 
two parts, i.e., symmetric eigenstates \s^), jss), \sq), jss), 
and antisymmetric eigenstates |s4) and jsy). By sym- 
metric we mean the state remains unchanged when one 
exchanges the particles 1 and 3, hence for symmetric 
state the entanglement between particles (1,2) and (2,3) 
equals. For the antisymmetric states, it is easy to check 
that these is no entanglement between particle (1,2) or 
(2,3). These observations for the eigenstates together 
with their distribution in the steady, we can conclude 
that the entanglement between the nearest neighbors are 
equal although the two temperatures for the two baths 
are different, this analysis confirms the numerical simu- 
lation presented in the figure. 



The top-right figure in Fig [5] shows the probability dis- 
tribution for the nonsymmetric case. We have set e=3, 
Ji=0.5 and J2=2.5. From this figure, we find that entan- 
glement in the steady states is mainly determined by the 
state \sq). We plot the probability distribution of \sq) 
with different mean bath temperature in the bottom of 
Figini Now we analyze the entanglement properties for 
this state. By tracing out one particle from the state |s6)j 
we can easily obtain the concurrence for the remaining 
two spins as Ci2=cos0, Ci3=cos0sin0 and C23=sin0. 
In nonsymmetric case, for example, Ji=0.5, J2=2.5, we 
have tan 0=5, which results in sin 6*^1. This is the reason 
why C12— Ci3 in Fig[4] and they are much smaller than 
C23. Moreover, observing the distribution probability for 
\sq) (blue line in the middle figure of FigOand the bot- 
tom figure), we can find that the peak of entanglement 
does not appear at Ar=0 in nonsymmetric case, namely 
the temperature difference favors the steady state entan- 
glement. 

Finally we study the effect of two coupling constant 
on the concurrence for fixed bath temperature. In FigjHl 
we show the concurrences as functions of Ji and J2 with 
£=3, rM=0.5 and AT= - 0.4. Both Ji and J2 enhance 
the entanglement. This enhancement is more strikingly 
in the case of Ji for C12 while J2 for C23. Due to the tem- 
perature difference, C13 is not symmetric about Ji—J2- 

In summary, we have studied the steady state entan- 
glement in a three-qubit XX model. The qubits are cou- 
pled to two independent bosonic baths at different tem- 
peratures. With the help of the effective Hamiltonian 
approach, we have calculated the steady state entangle- 
ment and discussed its dependence on temperatures and 
couplings. Two types of the coupling, i.e., symmetric and 
nonsymmetric one are considered. When the coupling is 
symmetric, we find that the entanglement between the 
nearest neighbors, i.e., C12 and C23, are equal though the 
temperatures of the two baths are different. The max- 
imum of entanglement is found when the temperatures 
of the two baths are equal. For the nonsymmetric case, 
however, the maximal entanglement arrives when the two 
baths at different temperatures. By analyzing the distri- 
bution for each eigenstate, we qualitatively explain these 
interesting phenomena. The dependence of the entangle- 
ment on the coupling constants are also presented and 
discussed. 
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